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We study a generally covariant model in which local 
Lorentz invariance is broken by a dynamical unit timelike vec- 
tor field It" — the "aether" . Such a model makes it possible to 
study the gravitational and cosmological consequences of pre- 
ferred frame effects, such as "variable speed of light" or high 
frequency dispersion, while preserving a generally covariant 
metric theory of gravity. In this paper we restrict attention 
to an action for an effective theory of the aether which in- 
volves only the antisymmetrized derivative V[aUb]- Without 
matter this theory is equivalent to a sector of the Einstein- 
Maxwell-charged dust system. The aether has two massless 
transverse excitations, and the solutions of the model include 
all vacuum solutions of general relativity (as well as other 
solutions). However, the aether generally develops gradient 
singularities which signal a breakdown of this effective the- 
ory. Including the symmetrized derivative in the action for 
the aether field may cure this problem. 



I. INTRODUCTION 

The Lorentz group is non-compact, since the boost pa- 
rameter is unbounded. This makes exact Lorentz invari- 
ance impossible to test uniformly. Lorentz invariance has 
thus been tested only up to some maximum boost and be- 
yond that lies an infinite volume of uncharted territory 
in the Lorentz group. Contrast this with the rotation 
group. Rotation invariance can be tested by filling in the 
compact 5*0 (3) group manifold more and more densely 
with data points, or by checking a few randomly selected 
rotations. The rotation group can be and has been uni- 
formly explored. 

There is also reason to doubt exact Lorentz invariance: 
it leads to divergences in quantum field theory associated 
with states of arbitrarily high energy and momentum. 
This problem can be cured with a short distance cutoff 
which, however, breaks Lorentz invariance. 

For these reasons we entertain the possibility that there 
is a preferred rest frame at each spacetime point. In par- 
ticular, we seek a viable effective field theory incorporat- 
ing a breaking of local Lorentz invariance. 

If the preferred frame were to be a fixed external struc- 
ture, then it would violate general covariance, which 
would require us to abandon general relativity (or any 



generally covariant modification thereof). General co- 
variance ordinarily implies that the divergence of the 
matter energy-momentum tensor Tab vanishes when the 
matter field satisfies its equation of motion. This is 
required for consistency of the Einstein field equation^ 
Gab — SnGTab, since the divergence of the Einstein ten- 
sor Gab is identically zero by virtue of the contracted 
Bianchi identity. If a fixed preferred frame is introduced 
into the matter action, for example, general covariance is 
lost. Tab is not divergenceless, and the Einstein equation 
is inconsistent. We therefore seek to incorporate the pre- 
ferred frame while preserving general covariance, which 
requires that the preferred frame be dynamical. 

It would be most appealing if the preferred frame were 
somehow determined by the spacetime metric itself. As 
discussed below, in spacetimes with an initial singularity 
the metric can be used to define a cosmological time func- 
tion, the gradient of which then determines a preferred 
frame (which is by construction timelike). However, the 
nonlocal relationship between this frame and the metric 
results in infinitely nonlocal field equations if this frame 
is incorporated into the action principle. Theories with 
even a finite amount of nonlocality are generally horribly 
unstable , so we do not consider this a viable approach. 

To avoid such unacceptable behavior the preferred 
frame should arise from local conditions, which of course 
reflect conditions at earlier times but only through dy- 
namics. For example, this dynamical frame could be 
defined by a vector field or by the gradient of 

a scalar field [^j7|. In these approaches the presence or 
absence of a preferred frame depends on the field con- 
figuration, since the preferred vector may vanish or may 
not be timelike. 

Since our motivation arises from doubts about the fun- 
damental validity of exact Lorentz invariance, we are 
not interested in a theory possessing a Lorentz invari- 
ant phase. We wish to study instead an effective theory 
in which there is always a preferred frame. This frame 
is defined by a timelike direction or, equivalently, by a 
unit timelike (contravariant or covariant) vector field.^ 
Such a field carries a nonlinear representation of the lo- 
cal Lorentz group since the field takes values not in a 
vector space but on the unit hyperboloid in the tangent 
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^We use units with c — 1 and the metric signature (H ). 

■^A unit timelike vector contains a discrete piece of informa- 
tion that a frame by itself does not have, namely, a time 
orientation. The theory may or may not depend on this 
orientation. 
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space. This could therefore be called a theory of non- 
linearly realized Lorentz invariance. It is analogous to a 
gauge theory with a non-linear sigma model Higgs field 
of fixed norm [^,1). There seems to be no generally ac- 
cepted terminology for this sort of symmetry structure. 
Since the symmetry breaking unit vector field is not a 
state-dependent expectation value but rather breaks this 
symmetry in all states, it may be misleading to say the 
symmetry is "spontaneously broken" . For lack of a bet- 
ter idea we shall take refuge in ambiguity and just call it 
"broken". What this really means is that in order to im- 
plement the local Lorentz symmetry one must transform 
not only the matter fields but also the background unit 
vector, so for practical purposes it will appear as what 
would normally be called broken Lorentz symmetry. 

The theory described here was devised about a decade 
ago by J.C. Dell together with one of us [0 , and we have 
since learned that similar ideas have independently been 
studied previously. The non-gravitational part of the the- 
ory (and generalizations thereof) was considered by Dirac 

in the early 1950's as a new theory of electrons (in 
which the unit timelike vector played the dual role of 
gauge-fixed vector potential and fiow vector of a stream 
of charged dust). A class of generally covariant theo- 
ries breaking Lorentz invariance was studied by Gasperini 

in many papers. In this work the tetrad formalism 
was used, and the local Lorentz symmetry was broken by 
including in the action terms referring to a fixed "inter- 
nal" unit timelike vector. This is equivalent to our formu- 
lation in terms of the metric and a unit timelike vector. 
To see the equivalence, note that the tetrad defines a 
metric and associates to the fixed internal vector a unit 
vector field on spacetime. The only other information in 
the tetrad is the gauge freedom parametrized by the local 
rotations leaving invariant the preferred timelike vector. 
Eliminating this gauge freedom leads to the formalism 
used in this paper. Gasperini has studied both cosmo- 
logical and central field solutions, with various choices 
for the specific form of the second derivative terms in the 
Lagrangian. In the present paper we focus on a different 
Lagrangian. 

The particular Lagrangian studied here was also con- 
sidered by Kostelecky and Samuel as a simplified 
model of the spontaneous Lorentz symmetry breaking 
that might occur in string theory, although in Ref. 
the preferred vector was not necessarily timelike. More 
generally, those authors argued that spontaneous Lorentz 
symmetry breaking in string theory may produce vacuum 
expectation values of more than one tensor field. In this 
case, rather than having just a single "preferred frame" 
there might be several background Lorentz tensors which 
collectively break part or all of the Lorentz symmetry. 

The remainder of this paper is organized as follows. 
In section ^ we explain the nonlocality problem encoun- 
tered if a cosmolo gica l time is used to define the preferred 
frame. In section 



ory comprise a subset of the solutions to the coupled 
Einstein-Maxwell-charged dust equations. Several types 
of exact solutions to the field equations are character- 
ized in section and the linearized theory is studied 
in section Coupling of the preferred frame to matter 
fields is discussed in section Both the dimension < 4 
couplings and some higher dimension ones are examined. 
The paper concludes with a brief discussion in section 



VII 



Among the dimension > 4 couplings are included the- 
ories involving Lorentz non-invariant dispersion at high 
wavevectors. This is motivated by recent work in which 
high frequency dispersion is invoked to avoid the role of 
trans-Planckian modes in the Hawking effect (for a re- 
view see [0). In this framework such theories can be 
formulated in a generally covariant manner so that grav- 
itational effects can be consistently incorporated. Higher 
dimension couplings also also provide an alternative gen- 
erally covariant formulation of variable speed of light 
models, in which different fields propagate at different 
speeds possibly at different cosmological epochs. Such 
models have recently been of interest as potential alter- 
natives to standard cosmology, and have been given gen- 
erally covariant formulations using additional fields to 
define the preferred frame 



II. COSMOLOGICAL TIMES 

In this section we briefly describe the construction of 
a cosmological time function determined purely by the 
metric, and the reason for rejecting it for the purposes of 
an effective theory of local Lorentz symmetry breaking. 

The cosmological metric of our universe, by virtue of 
its (approximate) homogeneity, (approximately) defines 
a preferred spacclike foliation of spacetime. However, 
this particular definition of the time function relies on 
the symmetry of the spacetime. For a workable theory 
with general covariance what is needed is a definition 
of cosmological time that can be used independently of 
symmetry. 

It is difficult to think of a notion of cosmological time 
that would make sense for all possible cosmologies. How- 
ever, if we restrict attention to spacetimes with a "begin- 
ning" , then two notions of cosmological time at x present 
themselves: (i) volume time, the spacetime volume (or 
perhaps the fourth root thereof) of the past set /"[x], 
and (ii) maximal time, the maximal proper time along a 
causal curve going back to the initial singularity.^ Other 



The maximal time function has been discussed in Refs. 



[II our proposed field theory of a pre- ^j^. 



ferred frame is formulated and its general properties are 
investigated. It is seen that the solutions for our the- 



[L5-|17|]. In particular, a powerful theorem proved in Ref. 
[ L7 establishes a number of properties of this function under 
the further assumption that the initial singularity is the only 
place past directed causal curves can end. 
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possibilities are combinations or smoothed averages of 
these times. Both of these time functions are determined 
non-locally but causally by the spacetime to the past of 
X. They may or may not be sufficiently smooth functions 
to enter meaningfully into a local action principle.^ If we 
assume that they are indeed sufficiently smooth we find 
that there is in any case a fatal problem with using them 
in this manner, as will now be explained. 

Suppose that to the usual action for gravity and matter 
fields is added a term involving one of the above cosmo- 
logical times, 



cal 



Sr. 



(2.1) 



We assume that the equations of motion are obtained 
as usual by requiring that the action is stationary with 
respect to variations of the fields. The variational deriva- 
tive of the action with respect to gab{x) has the form 



5S 



5S^ 



local 



lab{x) 5gab{x) 



6'Sr , / ,4 / 5t{x') 

+ / ax 



5t{x') 5gab{x) ' 



(2.2) 



where the prime on 5' indicates that the metric depen- 
dence of t[x) is not included in the variation. Since 
5t{x') / 5gab{x) has support when x' lies to the future of 
X, the field equation 5S / 5gab{x) — involves the val- 
ues of the fields to the future of x. Indeed the metric 
field equation is infinitely non-local in time, since the 
time function at any point to the future can be affected 
by a metric variation at x. Even finite nonlocality in 
time leads to unphysical instability [Q, so this approach 
to incorporating a preferred frame must be rejected. If 
the action depends on r only through its derivative VaT, 
then the equation of motion would be causal if (5Vct(x') 
depended only on gab{x) at a; = x' . However, this is not 
the case for the either the volume time or the maximal 
time. 



III. AETHER DYNAMICS 

We now turn to a class of theories in which there is a 
preferred frame which is determined by a local field. It 
is convenient to give a name to this field, and "aether" 
seems as good a name as any. Let us take the aether field 
to be a unit timelike vector field m", which is dimension- 
less, like the metric. To handle the condition that is 
a unit vector, we include in the action a Lagrange multi- 
plier term. Note that we are implicitly assuming that the 



*One of the results of the theorem of Ref. jT^ referred to 
in the previous footnote is that the maximal time function 
is locally Lipschitz and its first and second derivatives exist 
almost everywhere. The volume time function may well be 
even better behaved. 



spacetime admits a globally defined unit timelike vector 
field which is the case if and only if the spacetime is time 
orientable. 



A. Action 

The most general Lagrangian involving the metric and 
the aether with two or fewer derivatives is, up to a total 
divergence, 

Cg^u = ao - aiR - a2RabU°'u'' 

-bl F'^'Fab " 62 (VaUb)(V"«'') - 63 u'^ua, (3.1) 

where := u^Wmu"", and Fab is defined in analogy to 
the electromagnetic field strength. 



Fab := 2V[aU6]. 



(3.2) 



The term (VaU°)^ is equivalent, via integration by 
parts, to the combination (VaUfc)(V°M'') - (1/2) F"^ Fab + 



RabU^u'', so has not been included in (|3.lD. 

The Lagrangian ( ^T| ) is similar to the one discussed 
in Ref. ||l^] as the most general Lagrangian for a vector- 
tensor theory of gravity including terms up to second 
order in derivatives and quadratic in the vector field. The 
differences are that (i) the terms and Ru^ are missing 
from our action since the vector field is constrained to 
be a unit vector, and (ii) we have included the quartic 
term 6311^, which was omitted in Ref. because it is 
not quadratic in u". |l|l Note that, even without the last 
term, our theory is not a special case of the vector-tensor 
theories discussed in Ref. (l8| , since the constraint v? = 1 
affects the field equations. 



The coefficient ao in the action (3.1) has mass dimen- 
sion 4 while ai_2 and 61,2,3 have mass dimension 2. Lack- 
ing the underlying fundamental theory we do not try to 
assign a priori the values of these coefficients. A partial 
analysis of the observational consequences and limits on 
them has been done for the vector-tensor theories [p^ , 
however that analysis does not apply directly to our case 
due to the presence of the constraint term. It is fairly 
clear nevertheless that whatever values 61,2,3 take, agree- 
ment with observation will require that 02 ^ ai, and 
that ao/ai (which is basically the cosmological constant) 
must not be much larger than the squared Hubble con- 
stant. 

In this initial foray we shall restrict attention to the 
simple case in which the only terms with non-zero coef- 
ficients are R and F^ . That is, we set ao = 02 = 62 = 
63 = 0. The minimal theory we consider is thus defined 
by the action 

5'mm[gab,M", A] = j X ^/~^ 

-aiR - hi F'^^'Fab + \{gabu"-v!' - 1)) . (3.3) 
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This minimal theory is one of the models considered by 
Kostelecky and Samuel fl^ in the paper mentioned in 
section I. Those authors studied a broader class of mod- 
els in which is not necessarily constrained to be a 
unit vector but rather possesses a Lorentz-invariant po- 
tential energy with a minimum at some fixed value of 
Uau"' ■ They also allowed for extra, compact spatial di- 
mensions of spacetime, and examined cases where the 
symmetry breaking vector lies in the extra dimensions as 
well as cases where it lies in the four ordinary spacetime 
dimensions. Our paper by contrast is restricted to four 
dimensions and to a timelike vector of fixed norm. Later 
in this paper we shall also add matter terms to the ac- 
tion, including terms which couple the aether field to the 
matter. 

Note that Fab is invariant under the "gauge transfor- 
mation" 



Va/, 



(3.4) 



however the constraint v?' — 1 does not share this 
symmetry (nor do the additional couplings in general), 
so the theory is certainly not "gauge invariant". The 
constraint does have a limited version of this symme- 
try however, namely for those functions / satisfying 
{ua + VaDiu" + V/) = Uau" = 1. Thc general solution 
to this equation is 



where 



qab 



-gab + UaUb 



(3.5) 



(3.6) 



is the (positive definite) spatial metric orthogonal to 
Thus the action ( |3.3| ) is invariant under the gauge trans- 
formation ( |3.4| ) if / is chosen arbitrarily on a spacelike 
surface and then determined uniquely elsewhere (up to a 
discrete choice of sign) by integration of (3.5) along the 
flow of u"-. 



B. Field equations 

The equations of motion arising from the action ( |3.3| ) 
are 

Gab = - — iFaraFb"" " gab) + —UaUb, (3.7) 

ai ai 



ah 



2bi 



gabu^-u^ = 1. 



(3.8) 



(3.9) 



The metric equation (3/7) has the form of the Einstein 
equation Gab = ^T^GTab, where G = l/ldnai, and the 
stress tensor receives contributions from both thc F^ 



term and the co nstr aint term in the action. (The con- 



straint equation (3^) has been used to drop the contribu- 
tion to (3.7) that would have come from the variation of 
y/^g in the constraint term.) The contribution from the 
constraint term looks like that of a (pressurcless) dust 
with rest energy density 2A, and that from the F^ term 
is the usual Maxwell tensor familiar from electromag- 
netism, if we identify the vector potential as 



2\/biUr, 



(3.10) 



The stress tensor thus satisfies the usual energy condi- 
tions provided bi/ai and A/ai are positive 

In terms of th e ve ctor potential Am (3.10) the con- 
straint equation (3.9) becomes 



= 4&1 



(3.11) 



which can be interpreted as a gauge condition. The 
aether field equation ( |3.8| ) becomes the Maxwell equation 
with source equal to the current of a charged dust fluid 
with 4- velocity Ub and charge density —{X/^/bi). The 
evolution of A is determined by the current conservation 
equation which follows from divergence of the aether field 
equation ( |3.8[ ) upon using the identity VaVf.F'"'' = 0. 
Thus A satisfies a first order ordinary differential equa- 
tion along the flow lines of 



U^VaA = -AVait". 



(3.12) 



In particular, if A vanishes on a Cauchy surface, it must 
vanish everywhere. Also, thc sign of A on a given flow 
line cannot change, since if A = at any point on a flow 
line it must vanish everywhere on that line. 



Relation to Einstein-Maxwell-charged dust system 

We have just seen that the field equations of the mini- 
mal theory take the form of the coupled Einstein-Maxwell 
equations, with a charged dust matter source possessing 
charge to mass ratio — l/2\/5i. There is no explicit equa- 
tion of motion for the dust, however the normalization 

Taking the 



condition (3.9) provides such an equation, 
gradient of u'^ = 1 we have 

= Va{u\b) 
= 2u"S/aUb 
^2{u^VbUa+U^Fab). 



(3.13) 
(3.14) 
(3.15) 



Let us define Fab = 29(^^46] = 2^/hiFab- Then ( |3.15| ) 
becomes 



U VbUa 



1 



--Fr,bu\ 



(3.16) 



which is the equation of motion for a particle in the 
electromagnetic field Fab, with charge to mass ratio 
-l/2v^, the same ratio we inferred from the Einstein 
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and Maxwell equations! Thus any solution of our mini- 
mal theory is a solution of the Einstein-Maxwell-charged 
dust equations (although the converse is not true). The 
equivalence to a subset of the charged dust solutions 
demonstrates that the equations of our theory admit an 
initial value formulation, and it provides some useful in- 
tuition about the nature of the solutions. 

Our theory is not equivalent to the Einstein-Maxwell- 
charged dust system because in the general solution of 
that system the dust 4-velocity is not proportional to 
the vector potential in some gauge. That is, although 
there is always a gauge transformation that will make 
Am/'iy/bi a unit vector, it cannot in general be made to 
coincide with the dust 4-velocity.0 



IV. SOLUTIONS 

In this section we characterize a few types of solutions 
to the field equations. 

A. Solutions with A = 



If A = 0, then the two field equations (3.7,3 
the Einstein-Maxwell equations. Any solution to these 
equations is a solution in our theory provided a gau ge 
can be chosen so that the constraint equation ( [3.11 ) is 
satisfied. Such a gauge always exists, at least locally. 



Another way to think of this construction is in terms of 
the c ongru ence of integral curves of w". When Fab = 0, 
eqn. (3.15) implies that these curves are geodesies. More- 
over, if is the unit tangent field to a congruence of 
geodesies, then Fab — 2V[aUf,] = if and only if the 
congruence is hypersurface-orthogonal. Hence the gen- 
eral solution of this type is just an arbitrary solution to 
the Einstein equation, together with u° given by the unit 
tangent field of any hypersurface-orthogonal congruence 
of timelike geodesies in this metric. A special case is fiat 
spacetime, where the u"' congruence consists of straight 
lines normal to an arbitrary initial spacelike hypersur- 
face. 



Singular aether evolution 

This characterization of the Fab = solutions shows 
that, at least for such solutions, the evolution of Ua is gen- 
erally singular. The geodesies launched normally from a 
spacelike surface will typically cross. Where they do, the 
quantity VaUh will diverge. 

The existence of such singular evolutions for Ua signals 
a breakdown of the effective theory we are using. Perhaps 
are just it would be cured by including the term (V(QMf,))(V'-°u''-') 
in the action. (Without this term the action is insensi- 
tive to gradients for which the antisymmetrized deriva- 
tive V[aUb] vanishes.) For the purposes of the present 
paper we shall not pursue this question, but it should be 
addressed. 



B. Solutions with R, 







A special class of solutions to the field equations with 
A = are those with Fab = 0. For such fields, (3. S) 
implies that A = 0, and the field equations (3.7-3.9) re- 



duce to the ordinary vacuum Einstein equation together 
with the constraint — 1. When Fab = it follows, at 
least locally, that Ua — VaT for some function r, and the 
constraint then implies that VotV^t — 1. The general 
solution for such a function r can specified by assigning 
the value r = to an arbitrary spacelike surface, and 
determining r elsewhere by "uniform normal extension" , 
i.e. by the differential equation n"VaT = 1, where is 
the unit normal to the surface. 



^The general form of the discrepancy between these two 4- 
vectors was found by Dirac (see the second paper of Ref. ]ll|), 
who showed that (in four spacetime dimensions) there is al- 
ways a gauge in which Am./'2\^ = itm 4- fVm?), where ^ and 
rj are functions that are constant along the flow lines of w". 
Dirac included the functions ^ and rj as dynamical variables in 
order to obtain a theory in which arbitrary electron streams 
were admitted. In the third paper of Ref. he allowed for 
multiple streams. 



Cosmological solutions 

If u° shares the symmetry of a homogeneous isotropic 
cosmological metric, then Fab ~ 0. The presence of the 
aether field therefore has no influence on the cosmological 
evolution unless there are additional terms in the action 
beyond the minimal model. In Ref. |^0| we examine some 
cosmological effects of coupling to a scalar field tho ugh a 
fourth spatial derivative term as discussed in Sect. VI B. 



Black hole solutions 

For a spherically symmetric black hole, a suitable con- 
gruence of geodesies is given by the radial free-fall tra- 
jectories that all have the same Killing energy, i.e. the 
same asymptotic velocity at spatial infinity. The same 
construction can even be applied in the case of a Kerr 
black hole, at least for the geodesies that are at rest at 
spatial infinity. This follows from the work of Ref. , in 
which this congruence is employed to construct a coordi- 
nate system for the Kerr metric using the time function 
T mentioned above. 
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C. Spherically symmetric, static solutions 

Here we seek to characterize the general sphericahy 
symmetric, static solution. We shall find that, besides 
the mass, the metric in these solutions has an additional 
free parameter, the "aether charge" . 

Some of the linearized static, spherically symmetric so- 
lutions were previously studied in Ref. [|3|. Those au- 
thors examined the case where m° is spacelike, and while 
in four spacetime dimensions they restricted attention 
to vanishing Lagrange multiplier A and vanishing field 
strength Fab- In the present work we treat the nonlin- 
ear case, considering only timelike m° and imposing no 
further restrictions on the fields. 

Coordinates can be chosen so the line element takes 
the form 

ds^ = gttdt'^ + grrdr^ - r'^idO'^ + sm^edip"^), (4.1) 
and the aether field takes the form 



(O)), so that At = A't + l3 and Ar = A'^ + -f-r- The 



u — Ut{r)dt + Ur{r)dr. 



(4.2) 



The only potentially nonzero c omp onent of Fab is then 
Frt — drUt, and the constraint (|3.9|) implies 



rr 2 
9 "r 



1. 



(4.3) 



The aether field equation (3.8) in coordinate form 
reads 



1 



:9a {^g'-^g^'^F^,) = -^u\ (4.4) 



A 



or, taking into account the form of the metric (4.1) 



A 



(4.5) 



The left hand side vanishes when /? = r, hence the field 
equation implies that Ait^ — 0, which in turn implies 
that either = or A = 0. In the former case, u° is 
proportional to the timelike Killing field itself. There are 
thus two cases to consider. 

If A 7^ there are in fact no static solutions, unless 
the coefficients in the action are such that the charge 
to mass ratio of the dust is extremal. Recall that any 
solution to our theory is a solution to the charged dust 
theory. However under the influence of gravitational and 
electric forces, the non-extremal charged dust cannot re- 
main static, since there is no p ressure p^p^ . 

If A = then (cf. Sect. IV A ) these are just the spher- 
ically symmetric static Einstein-Maxwell solutions, i.e. 
the Reissner-Nordstrom solutions, in a spherically sym- 
metric, static gauge with fixed norm (3.11). Such a gauge 
always exists, at least locally. Consider the gauge trans- 
form Ay, — A'^ + of an arbitrary vector potential A'^. 
To maintain spherical symmetry and time independence 
we must have a = (3t -\- "f{r) (using the coordinates in 



normalization (3.11) then implies 



7,r 



= -A;± J5«[5«(AJ+/3)2-46i], 



(4.6) 



where we have used the fact that for the Reissner- 
Nordstrom metrics grr = —^/gtt- In any region where 
g** and At are bounded one can always choose /? large 
enough so that the radical is real, and then j{r) can be 
found by integration. 

If a horizon is present g** diverges and it is not clear 
from the preceding discussion whether the unit timelike 
gauge can be accessed in a smooth way across the hori- 
zon. Indeed it can, however the maximal extension of 
the region over which such a gauge can be accessed de- 
pends on the parameter bi. This can also be understood 
from the equivalence with a charged dust solution. The 
radial congruence must satisfy the Lorentz force equa- 
tion ( ^.16| ), and this congruence can be nonsingular and 
time-independent only if the trajectories are monotonic 
in the r coordinate. In general, however, the trajectories 
bounce inside the black hole. 



Comparison with observation 

To compare with observation it would be necessary 
to determine which of the above solutions to use in the 
presence of a spherically symmetric static source such as 
a planet, star, or black hole. The metric associated with 
one of these objects depends on its "aether charge" — 
the charge of the Reissner-Nordstrom solution — which is 
determined by the "charge" of the "aether dust" that fell 
in when the object condensed. The choice is determined 
by the initial conditions on A, which are presumably cos- 
mological in origin. We have no theory of these initial 
conditions at this stage, but agreement with observations 
can put a bound on the amount of aether charge. If this 
charge is zero, then we have the usual Schwarzschild so- 
lution of general relativity (and, as discussed in section 
[VB, the aether field is the tangent field to a hypersur- 
face orthogonal congruence of timelike geodesies), which 
of course agrees with observations. 



V. LINEARIZED THEORY 

In this section we study the linearized equations de- 
fined by expanding about a background solution. 



(0) I u 
gab = gab + hab 



Ua = 



A = A(o)+A(i). 



For the background we take the flat metric 
,(") 



(0) 



(5.1) 
(5.2) 
(5.3) 

= rjab 



and a constant Ua ■ In this background solution the 
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equations of motion imply that the Lagrange multipher 
A*-°^ must vanish, hence we shall use the letter A for the 
perturbation A'"'^'. In this section we use the flat back- 
ground metric to raise and lower indices. Note that we 
use the perturbation of the covariant vector Ua to define 
the perturbation Va- 

We choose cartesian coordinates (a;°, a;*), i ~ 1, 2, 3, in 
which the components of rjab are diag(l, —1, —1, —1) and 

those of ui"' are (1,0,0,0). The linearized field equa- 
tions for this theory were first written down in a gen- 
eral gauge in Ref. who also pointed out that the 
Lorentz gauge can be accessed using the linearized dif- 
feomorphism invariance of the action. The Lorentz gauge 
condition for the metric perturbation is hab'^ = 0, where 
hab = hab — \hriab is the trace- reversed metric perturba- 
tion. In this gauge, the linearized equations of motion 
are 



^ (0) (0) 



2—u 

ai 



aVb-db{d''Va)=-^'> 
Zui 



,(0) 



habU 



vb = 



(5.4) 



(5.5) 



(5.6) 



In a source-free region, the residual gauge freedom is 
usually employed to set /iqi — h^i — and h = —h = Q. 

The possibility of doing so depends on the fact that 
these quantities satisfy the wave equation in general rel- 
ativity. In our case, the components /loi satisfy the wave 
equation, however h do es n ot, due to the source term on 
the right hand side of (5.4) that is always present (even 
outside matter) unless A = 0. This source corresponds 
to the energy density of the "charged dust" , and we wish 
to allow for the presence of this term. Therefore, rather 
than setting h to zero, we choose to set to zero the trace of 
the spatial part, h\, which does satisfy the wave equation. 
The proof that this can be done follows the same logic as 
in the usual case. This gauge condition implies h = hoo, 
hence /iqo — ^oo — 5^*700 — ^^qq, and the Lorentz gauge 
condition implies 9o^oo = and dihij = 0. Thus hoo is 
time-independent and the spatial part hij is a transverse 
traceless solution to the wave equation. This is not quite 
the same as the usual transverse traceless gauge in gen- 
eral relativity however, since hij = hij — ho^Sij, so hij is 
not transverse unless dihoo — 0, and h = — 2/ioo 7^ 0. 

It remains to consider the linearized equations for hoo, 
Vg, and A. The 00-component of the metric equation 
(5.4) and the constraint equation ( ^.6[ ) determine A and 
vo in terms of hoo- 

A -aiV^/ioo (5.7) 
Vo = -5/100- (5.8) 



Using the time- independence of vo , the time and space 
components of the aether equation (|5.5|) read 



V\o-5o(a'«,)--7^ 
Dv,~d,id^v,)^0. 



(5.9) 
(5.10) 



Let us decompose Vi into transverse and longitudinal 
parts, Vi = vf + vf , where d''vj — and vf = dif for 
some scalar field /. Then ( ^.10| ) implies that the trans- 
verse part satisfies the wave equation, Ovf = 0, so the 
the aether field has two transverse massless modes. 
As for the longitudinal part , equation (5^) implies 



doid'v^ 



bi - Qi 
2bi 



V^h 



00, 



(5.11) 



so in particular vf" has at most linear time dependence. 
(The same conclusion follows from the divergence of Eq. 
(P^.) Thus uf d,a{x)t + dib{x), where a{x) is de- 
termined by hoo Sind b{x) is arbitrary. 

In summary, the perturbation spectrum consists of two 
massless transverse traceless modes of hij , and two mass- 
less transverse modes of Vi. In addition there is a mode 
in which hoo is an arbitrarily specified time-independent 
function, which determines uo, A, and the time derivative 
of the longitudinal part wf. The time-independent part 
of Vi is also arbitrary. This last freedom corresponds to 
the linearization of the restricted gauge symmetry ( ^3.51 ). 

The longitudinal mode looks very strange at first sight. 
In the charged dust interpretation, the dust energy den- 
sity is adjusted to produce an arbitrary gravitational po- 
tential hoo, smd the perturbed metric, electromagnetic 
field, and charge density are all time independent, while 
the perturbed dust world lines are time dependent. This 
is a peculiarity of the first order perturbative solution 
however. No exact solution shares this property, as can 
be easily seen from the aether field equation (3.8). If 



the left hand side is invariant with respect to a time- 
like Killing field, and if A is also invariant, then so must 
be u''. Evidently the higher order terms in the equa- 
tions of motion induce time dependence into the solution. 
A similar phenomenon can be seen upon expanding the 
simpler Einstein-neutral dust system about the flat space 
solution with constant dust 4- velocity and vanishing den- 
sity. As in our case, the dust density perturbation can 
set up any static metric perturbation, and the linearized 
geodesic equation for the dust yields a time-dependent 
dust velocity perturbation. 



VI. MATTER COUPLINGS 

We have so far considered only the terms in the action 
involving the metric and the aether field and up to two 
derivatives. Suppose a matter term Smat[gab,u°',''P] is 
adde d to the action Smin [dab, u"", A] of the minimal theory 
( |3.3D , where -0 stands for a generic matter field. The 
variation of Smat with respect to the metric produces an 
additional contribution to the stress-energy tensor, and 
the variation with respect to u° produces an additional 
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term in the current on the right hand side of (3.8). The 
resulting field equation takes the form 



The identity VV''Fab = then impHes 

2 Su" 



(6.1) 



(6.2) 



which shows that now even if A is initially zero it need 
not remain zero. In the presence of such matter couplings 
the equivalence to charged dust is lost. 

We now consider specific types of matter couplings, 
first of dimension less than or equal to four, and next of 
dimension greater than four. 



The indices /, J in ( |6.3D are generational indices, and 
the coupling constants {clUj, cb, and c$ are all dimen- 
sionless. The ellipses in (S.3) stan d for similar terms for 
the other fermions, while those in (6^) stand for similar 
terms for the other gauge fields. 

Such additional kinetic terms modify the propagation 
speed of the various fields. For example, the propagation 
speed for the Higgs field, with respect to the preferred 
frame, is (l+Ccj.)"^/^, which is less or more than the speed 
of light if c$ is positive or negative respectively. The 
coupling constants must therefore be small numbers for 
fields whose propagator has been measured accurately. 
It would be interesting to determine what limits can be 
placed on these coefficients, particularly for fields such as 
the Higgs or gluons whose propagators are presumably 
not yet so well measured. 



A. Couplings of dimension < 4 



B. Couplings of dimension > 4 



A complete classification of Lorentz violating, gauge- 
invariant extensions of the SU(3)xSU(2)xU(l) minimal 
standard model has been given by CoUaday and Kost- 
elecky p5| , restricting attention to operators whose mass 
dimension is less than or equal to four, so as to preserve 
power-counting renormalizability. This class of low en- 
ergy effective actions includes both C PT-even and CPT- 
odd terms, and involves various coupling tensors with 
"generation" indices allowing for mixing of fermions from 
different generations. These coupling tensors are sup- 
posed to be Lorentz violating vacuum expectation val- 
ues arising in a theory with a fundamental underlying 
Lorentz symmetry. 

Here we consider the above class of Lorentz violating 
terms, keeping only those couplings that can be con- 
structed with the aether field u". With this restriction 
the antisymmetric tensor couplings are excluded, which 
rules out Lorentz-violating Yukawa couplings and cou- 
plings of gauge field strengths to Higgs bilinears, and 
limits the form of modifications of the gauge field kinetic 
terms. Invariance under time reversal w° — would 
be required if the physical significance of the aether is 
only to define a preferred frame and not a preferred lo- 
cal time orientation. If we accordingly further assume 
this symmetry, all the CPT-odd terms are excluded, 
which rules out terms with a vector coupled to fermion or 
Higgs currents, gauge field Chern-Simons currents, and 
the U(l) potential. 

The only possibilities remaining after all these restric- 
tions have been imposed are the modifications of the 
fermion, gauge field, and Higgs kinetic terms: 



\i{cL)i,ju''u''LnaDbLj + h.c. 



jCbU U g rlam^t 



bn 



(6.3) 



(6.4) 



(6.5) 



Once the restriction to terms of dimension 4 or less is 
dropped, the possibilities for Lorentz violating terms — 
like those for Lorentz invariant ones — are endless. Here 
we would like to consider just two types, which illustrate 
different possibilities that arise in the presence of Lorentz 
symmetry breaking. 



Modified kinetic terms 

If the couplin g coeffic ient for a Lorentz violating ki- 
netic term like (3.3-S.5) is field dependent and polyno- 
mial, rather than a constant, then the term is a dimen- 
sion > 4 operator. In this case it is possible that the 
coefficient was larger in the early universe than it is to- 
day, due to the cosmological evolution of the field(s) on 
which the coupling function depends. This provides an 
alternate approach to constructing generally covariant, 
variable speed of light cosmologies. Approaches using a 
vector 1^ or a scalar to define the preferred frame 
have been the subject of some recent papers. 



Modified dispersion 

Next we consider a deviation from Lorentz invariance 
that becomes strong only at high wavevectors. In the 
early universe, when the fields were highly excited at 
large wavevectors, the gravitational effects of such a de- 
viation could have been of paramount importance. The 
study of a model incorporating such effects is left to an- 
other paper [ po| . Here we indicate only an example of 
a term in the Lagrangian that produces high frequency 
dispersion in the propagation of a matter field, and we 
display the the form of the resulting contribution to the 
energy-momentum tensor. 
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Consequences of non-Lorentz invariant high frequency 
dispersion for the Hawking effect have previously been 
studied using 1+1 dimensional model field theories in 
which higher spatial derivative terms are added to the 
action (for a review see |l^), and recently such models 
have been generalized to field theory in the background of 
a 3+1 dimensional Robertson- Walker spacetime in order 
to study the consequences for the spectrum of primordial 
density fluctuations in inflationary cosmology . These 
models can be extended to an arbitrary 3+1 dimensional 
setting, preserving general covariance as well as spatial 
rotation symmetry in the local preferred frame. As an 
example consider the Lagrangian 



£^-i(vVVa^ + fco-'pV)')- 



(6.6) 



Here fco is a constant with the dimensions of inverse 
length which sets the scale for deviations from Lorentz 
invariance, and is the covariant spatial Laplacian, i.e., 



(6.7) 



where Da is the spatial covari ant derivative operator ||24 
and qab is the spatial metric (^^). 



The w°-dependence of the Lagrangian (6.6) produces 
a "matter" term in the aether field equation (6J_). The 
energy-momentum tensor for this Lagrangian is 

Tab = Va</?Vf,(p - gab 

2D\u"'u(^ay\m\Db)f + 2V„iiD^^q(^a''')yb)^ 

-V"\qabD^ipDmip)\. (6.8) 



In Ref. (20| we evaluate the expectation value of this 
energy-momentum tensor in a thermal state in flat space- 
time, which allows us to determine the modification of the 
equation of state produced by the fourth derivative term. 
This equation of state is then be used to study how the 
cosmological evolution is affected by the high frequency 
dispersion. 



VII. DISCUSSION 

We have made an initial attempt to study the possi- 
ble consequences of incorporating a preferred frame — the 
aether — into a generally covariant theory. With the ac- 
tion adopted in this paper the aether vector generically 
develops gradient singularities even when the metric is 
perfectly regular. We take this as a sign that the the- 
ory is unphysical as an effective theory (although if the 
aether sector is ignored the theory can be made to agree 
with observations with an appropriate choice of initial 
conditions, i.e. by setting Fab to zero.) The primary 
open questions are (i) what determines the initial values 
of the aether field and the Lagrange multiplier field, and 
(m) are the gradient singularities, which appear to be 



generic in the evolution of the aether, eliminated by in- 
cluding a symmetrized derivative term (V(aWf,))(V^''u''') 
in the action along with the antisymmetrized derivative 
term used in this paper? It is plausible that adding the 
symmetrized derivative term will have a significant effect, 
since with it the action is sensitive to the existence of any 
large gradients. 
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